NAVAL  POSTGRADUATE  SCHOOL 

Monterey ,  California 


THESIS 


ENHANCEMENT  OF  BOILING  HEAT  TRANSFER 
IN  DI-ELECTRIC  FLUIDS 

by 

Robert  A.  Egger 
September  1991 


Thesis  Advisor 


M.D.  Kelleher 


Approved  for  public  release;  distribution  is  unlimited 


T259716 


;lassified 

nty  classification  of  this  page 


REPORT  DOCUMENTATION   PAGE 


Report  Security  Classification  Unclassified 


lb  Restrictive  Markings 


security  Classification  Authority 


Declassification/Downgrading  Schedule 


3  Distribution/Availability  of  Report 

Approved  for  public  release;  distribution  is  unlimited. 


erforming  Organization  Report  Number(s) 


5  Monitoring  Organization  Report  Number(s) 


Name  of  Performing  Organization 
ival  Postgraduate  School 


6b  Office  Symbol 
(if  applicable)  34 


7a  Name  of  Monitoring  Organization 
Naval  Postgraduate  School 


Address  (city,  state,  and  ZIP  code) 
onterey,  CA  93943-5000 


7b  Address  (city,  state,  and  ZIP  code) 
Monterey,  CA  93943-5000 


Name  of  Funding/Sponsoring  Organization 


8b  Office  Symbol 
(if  applicable) 


9  Procurement  Instrument  Identification  Number 


Address  (city,  state,  and  ZIP  code) 


10  Source  of  Funding  Numbers 


Program  Element  No    Project  No     Task  No     Work  Unit  Accession  No 


Title  (include  security  classification)  ENHANCEMENT  OF  BOILING  HEAT  TRANSFER  IN  DI-ELECTRIC  FLUIDS 


Personal  Author(s)  Robert  A.  Egger 


l  Type  of  Report 
aster's  Thesis 


13b  Time  Covered 
From  To 


14  Date  of  Report  (year,  month,  day) 

September  1991      


15  Page  Count 
94 


Supplementary  Notation  The  views  expressed  in  this  thesis  are  those  of  the  author  and  do  not  reflect  the  official  policy  or  po- 
Ion  of  the  Department  of  Defense  or  the  U.S.  Government. 


Cosati  Codes 


Id 


Group 


Subgroup 


18  Subject  Terms  (continue  on  reverse  If  necessary  and  Identify  by  block  number) 
word  processing,  Script,  GML,  text  processing. 


Abstract  (continue  on  reverse  if  necessary  and  identify  by  block  number) 

Direct  application  of  two-phase  heat  transfer  in  the  liquid  cooling  of  electronic  components  in  flouorinated  hydrocarbons 
C-72),  is  severely  inhibited  by  the  excessive  amount  of  superheat  required  to  initiate  nucleate  boiling.  To  conduct  an  ex- 
rimental  study  of  nucleate  pool  boiling  of  FC-72,  an  experimental  test  chamber  was  constructed.  This  chamber  utilized  five 
rizontal  platinum  wires  of  0.05  mm  diameter  spaced  2.0  cm  vertically  from  each  other.  The  lowest  wire  was  progressively 
ated  from  the  natural  convection  region  through  nucleate  boiling,  and  a  study  was  made  on  the  effects  of  the  boiling  wake 
ime  on  the  heat  transfer  rate  of  the  upper  wires. 


Distribution/Availability  of  Abstract 
unclassified/unlimited  D  same  as  report 


□  DTIC  users 


21  Abstract  Security  Classification 
Unclassified 


i  Name  of  Responsible  Individual 
•D.  KeUeher 


22b  Telephone  (include  Area  code) 
(408)  646-2530 


22c  Office  Symbol 

ME/Kk 


'  FORM  1473,84  MAR 


83  APR  edition  may  be  used  until  exhausted 
All  other  editions  are  obsolete 


security  classification  of  this  page 
Unclassified 


Approved  for  public  release;  distribution  is  unlimited. 

Enhancement  of  Boiling  Heat  Transfer 
in  Di-Electric  Fluids 


by 


Robert  A.  Egger 

Lieutenant  ,  Unitea  States  Navy 

B.E.  Mechanical  Engineering,  Cleveland  State  University,  1985 

Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degrees  of 

MASTER  OF  SCIENCE  IN  MECHANICAL  ENGINEERING 

and 
MECHANICAL  ENGINEER 

from  the 

NAVAL  POSTGRADUATE  SCHOOL 
September  1991 


ABSTRACT 

Direct  application  of  two-phase  heat  transfer  in  the  liquid  cooling  of  electronic 
components  in  flouorinated  hydrocarbons  (FC-72),  is  severely  inhibited  by  the  excessive 
amount  of  superheat  required  to  initiate  nucleate  boiling.  To  conduct  an  experimental 
study  of  nucleate  pool  boiling  of  FC-72,  an  experimental  test  chamber  was  constructed. 
This  chamber  utilized  five  horizontal  platinum  wires  of  0.05  mm  diameter  spaced  2.0 
cm  vertically  from  each  other.  The  lowest  wire  was  progressively  heated  from  the  na- 
tural convection  region  through  nucleate  boiling,  and  a  study  was  made  on  the  effects 
of  the  boiling  wake  plume  on  the  heat  transfer  rate  of  the  upper  wires. 
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NOMENCLATURE 


Symbols 

Units 

Description 

A 
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Surface  Area  of  Test  Surface 
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Thermal  Expansion  Coefficient 
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J/kg 
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m/sec2 
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h 
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k 
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Thermal  Conductivity 

L 

m 
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T 

°K 

Temperature 

Ra 
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ND 

Prandlt  Number 

Nu 
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Nusselt  Number 
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sub 
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liquid 

vapor 

Saturation  Conditions 

Subcooled  Conditions 

Bulk  Fluid  Conditions 

Platinum  Wire 

2C1  Resistor 

Surface 


GREEK  SYMBOLS 
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Ns/m2 

Dynamic  Viscosity 

V 

m2  /s 

Kinematic  Viscosity 

o 

N/m 

Surface  Tension 

P 

kg/m2 

Fluid  Density 

CO 

Uncertainity  in  measurements 

I.     INTRODUCTION 

A.  HISTORICAL 

Ever  since  the  development  of  the  first  practical  solid  state  integrated  circuit  in  1959, 
there  has  been  a  continuous  emphasis  on  packaging  an  ever  higher  density  of  logic  cir- 
cuits on  a  single  substrate.  Corresponding  with  this  increase  packaging  density  is  an 
increasing  power  dissipation  requirement  along  with  an  increase  difficulty  in  the  thermal 
management  capabilities.  Initially,  thermal  management  techniques  consisted  of  noth- 
ing more  than  a  heat  sink  and  natural  convection  air  cooling.  This  naturally  progressed 
into  heat  sinks  with  forced  convection  air  cooling,  and  with  the  construction  of  large 
mainframe  computers  and  first  generation  super-computers  the  use  of  indirect  liquid 
cooling  was  initiated.  The  next  practical  method  for  heat  removal  is  direct  liquid  cooling 
or  immersion  cooling  of  microelectronic  components  using  dielectric  liquids. 

Although  direct  liquid  cooling  has  been  successfully  employed  for  nearly  40  years  in 
large  hi-powered  microwave  equipment,  klystron  tubes,  travelling-wave  guide  tubes,  and 
power  transistors  [Ref.  1  :  p.  12],  only  recently  has  it  been  employed  for  use  in  the  ther- 
mal management  of  microelectronics.  The  first  computer  systems  to  successfully  employ 
immersion  cooling  of  microelectronics  were  the  CRAY-2  and  the  ETA-10  supercom- 
puter [Ref.  2  :  p.  8].  Today  the  research  in  direct  immersion  cooling  is  continuing  with 
particular  emphasis  on  utilizing  dielectric  liquids  as  a  cooling  medium,  in  particular  the 
preffuorinated  FC  liquid  family  of  the  3M  corporation  [Ref.  3]. 

B.  BACKGROUND 

The  advantages  offered  by  direct  immersion  cooling  is  self  evident  from  Figure  1 
By  utilizing  forced  convection  direct  liquid  cooling  with  FC-77,  the  CRAY-2  has  been 
able  to  reduce  it's  required  volume  by  an  order  of  magnitude  as  compared  to  it's  indirect 
liquid  cooled  competitors.  This  was  being  achieved  in  the  CRAY-2  while  at  the  same 
time  maintaining  the  same  power  dissipation  of  the  competition,  [Ref.  1]. 
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Figure  1.      Power  consumption  vs  Volume:     Of  present  computers 


As  previously  stated,  the  CRAY-2's  packaging  density  was  achievable  due  to  direct 
liquid  immersion  cooling,  however,  it  did  use  a  forced  convection  cooling  scheme.  A 
thermal  management  scheme  that  relies  on  boiling  heat  transfer  would  provide  an  even 
more  efficient  heat  transfer  capability,  and  has  been  the  attention  of  considerable  re- 
search in  recent  years.  Unfortunately,  these  prefluorinated  dielectric  fluids  possesses 
anomalies  which  have  held  back  their  widespread  use  and  acceptance  of  direct  immersion 
cooling. 

In  comparison  to  water  and  many  other  coolants,  these  prefluorinated  dielectric 
liquids  posses  a  lower  value  for  their  thermal  conductivity,  and  heat  of  vaporization, 
they  also  have  an  extremely  small  surface  tension  which  translates  into  a  near  zero  wet- 
ting angle  and  therefore,  a  high  wettability  on  most  surfaces.  This  last  feature  results  in 
the  elimination  of  many  potential  nucleation  sites  which  in  turn  results  in  a  incipience 
super  heat  excursion  which  results  in  a  very  high  wall  temperature. 

In  this  study,  heat  transfer  enhancement,  and  therefore  a  surface  temperature  re- 
duction, due  to  bubble  pumping  in  FC-72  will  be  investigated. 

C.     PREVIOUS  WORK 

An  extensive  amount  of  research  has  been  conducted  in  the  boiling  regime  of 
dielectric  liquids.   The  most  helpful  of  the  literature  used  will  be  noted. 


You,  Simon  and  Bar-Cohen  [Ref.  4  and  Ref.  5  ]  conducted  an  analyses  on  the  ef- 
fects of  pressure,  subcooling,  and  dissolved  gas  content  on  the  boiling  heat  transfer,  and 
the  boiling  incipience  of  dielectric  fluids.  In  their  experiment  they  used  FC-72  as  the 
working  fluid  and  a  0.1,um  thin  film  of  platinum  on  a  0.51mm  diameter  quartz  cylinder, 
with  the  platinum  film  acting  as  both  heater  and  temperature  sensor.  They  concluded 
that: 

1.  An  increase  in  pressure  of  56kPa  reduces  wall  superheat  by  approximately  3°C  in 
the  fully  developed  nucleate  boiling  region. 

2.  Subcooling  has  little  effect  on  the  fully  developed  nucleate  boiling  region. 

3.  Nucleate  boiling  hysteresis  is  observed  for  both  the  subcooled  and  gas-saturated 
cases.   The  hysteresis  decreases  as  the  gas  content  increases. 

4.  Incipience  bubble  radius  computed  for  FC-72  was  0.05-0. 1/rni. 

5.  Effects  of  dissolved  gases  in  the  fluid  is  seen  at  low  content  values  as  an  increase 
in  run  to  run  variability. 

6.  The  effect  of  subcooling  on  pool  boiling  incipience  is  apparently  small. 

7.  A  wide  variation  in  incipience  superheat  values  for  nominally  identical  runs  were 
observed. 

You,  Simon,  Bar-Cohen  and  Tong  [Ref.  6]  conducted  an  experimental  investigation 
of  boiling  in  R-113.  The  experiment  was  conducted  using  0.13mm  diameter  Chromel 
wires  and  0.51mm  diameter  platinum  thin  film  heater.   The  following  was  found. 

1.  A  Scanning  Electron  Microscope  (SEM)  analysis  was  conducted  on  the  surface  of 
the  thin  film  heater,  in  which  the  chromel  wire  surface  was  found  to  contain  a 
higher  density  of  nucleation  sites  compared  to  the  thin  film  platinum  heater. 

2.  The  platinum  heater  surface  was  found  to  have  a  higher  incipience  superheat  than 
the  chromel  wires  ,  73°  C  versus  38°  C  maximum.  This  tends  to  correlate  well  with 
item  1. 

3.  Again,  a  wide  variation  in  incipience  superheat  values  for  nominally  identical  runs 
were  experience. 

Marto  and  Lepere  [Ref.  7]  conducted  experiments  with  surface  enhancements  for 
improving  pool  boiling  performance.  Marto  and  Lepere  using  FC-72  and  R-113  as 
working  fluids  compared  the  performance  of  three  commercially  available  boiling  en- 
hancement surfaces  to  plain  copper  tubing.  In  each  case  the  enhanced  surface  resulted 
in  a  lower  wall  surface  temperature  compared  to  the  plain  copper  tube.  However,  the 
required  incipience  superheat  for  nucleation  remained  unaltered. 

Bergles  and  Kim  [Ref.  2]  analyzed  methods  to  reduce  the  temperature  overshoot  of 
immersion  cooling.    Their  method  used  an  active  bubble  generator  in  the  fluid  located 


below  the  simulated  micro  chip.  They  used  the  bubble  generator  to  activate  dormant 
and  extinct  nucleation  sites.  They  obtained  a  reduction  of  incipience  superheat  required 
from  22-33°  C  to  8-12°  C  in  a  working  fluid  of  R-113. 

Maddox  and  Mudawar  [Ref.  8]  studied  the  effect  of  subcooling  and  surface  aug- 
mentation on  values  of  Critical  Heat  Flux  (CHF).  They  found  that  the  CHF  of  FC-72 
can  achieve  a  level  greater  then  360  Watts/cm2  with  combined  use  of  subcooling  and 
surface  augmentation. 

Danielson,  Tousignat  and  Bar-Cohen  [Ref.  9]  conducted  an  analysis  of  a  large 
number  of  commercially  available  prefluorinated  inert  fluids.  The  method  of  investi- 
gation used  a  0.25mm  diameter  platinum  wire  heating  element,  with  the  working  fluid 
at  saturation  conditions.  Important  points  reported  are: 

1.  The  critical  heat  flux  of  FC-72  was  reported  as  approximately  20.3  Watts/cm2. 

2.  The  nucleate  boiling  characteristic  of  highly  wetting  fluids  can  be  predicted  using 
correlations  developed  for  water  and  other  fluids,  but  are  highly  dependent  of  the 
heating  surface  condition. 

3.  Single  phase  convection  data  are  predicted  accurately  by  a  correlation  developed 
by  Kuehn  and  Goldstein.[Ref.  10] 

4.  The  superheat  excursion  at  incipience  boiling  always  appears,  but  the  magnitude 
varies  widely. 

5.  The  apparent  radii  of  nucleation  bubbles  were  calculated  to  be  0.1  to  0.2^m. 

Kuehn  and  Goldstein  [Ref.  10]  has  obtained  a  correlation  of  natural  convective  heat 
transfer  for  infinite  cylinders.  This  correlation  was  used  in  this  thesis  data  analysis  due 
to  the  low  rayleigh  numbers  obtained  in  the  experimental  runs. 

From  the  3M  Product  Manual  [Ref.  3]  it  can  be  found  that  the  highly  wetting 
dielectric  fluids  have  a  high  solubility  for  gases.  FC-72  can  contain  up  to  48%  air  by 
volume  compared  to  1.9%  air  for  water.  This  is  the  necessity  of  thoughly  degassing  the 
fluid  prior  to  the  experimental  runs. 

D.     OBJECTIVES 

The  objective  of  this  work  was  to  investigate  the  effects  of  the  boiling  wake  from  a 
heated  wire  upon  the  heat  transfer  rate  from  wires  placed  within  the  wake.  This  work 
was  carried  out  using  the  dielectric  fluid  FC-72. 


II.     EXPERIMENTAL  APPARATUS 

A.     DESCRIPTION  OF  COMPONENTS 

The  experimental  apparatus  used  in  this  study  is  a  redesign  of  an  experimental  ap- 
paratus used  in  a  previous  thesis.  Emphasize  on  this  design  was  to  minimize  systematic 
errors,  minimize  uncertainty,  and  to  correct  deficiencies  to  ensure  an  overall  increase  in 
the  data's  creditability  over  the  previous  apparatus,  [Ref.  11]. 

In  Figure  2  an  overall  view  of  the  experimental  setup  is  shown,  identifying  all  major 
components  that  were  utilized.  Figure  3  is  a  photograph  of  the  actual  experimental 
chamber. 

The  apparatus  chamber  consists  of  a  box  constructed  of  0.75  inch  thick, 
polycarbonate  sheets  with  a  0.25  inch  thick  aluminum  coverplate  fitted  with  a  rubber 
O-Ring  seal.  Inside  the  box  are  two  bulk  heaters  to  raise  the  fluid  to  saturation  tem- 
perature for  degassing  purposes,  and  to  maintain  the  bulk  fluid  temperature  near  con- 
stant. The  aluminum  cover  plate  is  mounted  with  three  thermo-electric  coolers,  and 
serves  as  a  condenser,  where  by  the  vapor  condenses  and  returns  to  the  chamber  by 
gravity.  An  insert  board  is  held  in  the  chamber  on  which  are  mounted  the  test  wires, 
tension  springs,  thermocouples,  and  voltage  pickofTs. 
1.     Test  Chamber 

The  interior  dimensions  of  the  chamber  are  6.0  inches  wide,  by  2.375  deep,  by 
6.25  inches  high  ,  as  shown  in  the  isometric  view  of  Figure  8.  The  entire  chamber  was 
assembled  using  IPN  #16  Acrylic  bonding  adhesive  and  a  bead  of  Di-Metheyl  Ethelene 
was  used  to  seal  all  seams  along  the  interior  surface  to  ensure  liquid  tightness. 

To  permit  access  for  instrumentation  and  power  leads,  and  as  a  means  to  ensure 
pressure  equalization,  a  0.75  inch  diameter  and  0.25  inch  diameter  hole  were  placed  into 
the  rear  face  of  the  chamber  at  the  lower  right  and  upper  right  corner  respectively,  as 
shown  in  Figure  4.  Mounted  securely  around  these  access  holes  are  lengths  of  tygon 
tubing  of  approximately  18  inches  in  length,  with  the  opened  ends  suspended  well  above 
the  chamber. 

On  the  interior  side  walls,  slots  are  milled  to  permit  the  placement  and  removal 
of  an  instrumented  wire  board  on  which  the  experiments  are  conducted.  Figure  6  is  a 
photograph  of  the  insert  board,  while  a  schematic  drawing  is  shown  in  Figure  7. 


2.  Aluminum  Cover  Plate 

A  0.25  inch  thick  by  7.00  inch  long  by  3.25  inch  wide  aluminum  plate  was  uti- 
lized for  use  as  a  cover  plate  for  the  chamber,  and  as  a  condenser  surface  for  the  FC-72 
vapor.  Mounted  on  the  interior  face  are  two  thermocouples  to  monitor  the  condensor 
surface  temperature,  as  shown  in  Figure  9.  On  the  exterior  face  of  the  plate,  three 
thermo  electric  coolers  were  installed  for  heat  removal.  To  ensure  a  proper  seal  between 
the  aluminum  plate  and  the  box  structure  ,  a  shallow  O-Ring  groove  was  machined  into 
the  interior  face  and  an  O-Ring  of  0.0624  inch  diameter  was  installed.  To  adequately 
secure  the  cover  plate  to  the  chamber  box,  eight  stainless  steel  screws  of  1/16  inch  di- 
ameter are  used  to  compress  the  O-Ring  against  its  mating  surface. 

As  mentioned  previously,  three  solid  state  thermoelectric  cooling  devices  from 
the  MELCOR  corp  were  employed  to  cool  the  cover  plate  condensor.  They  were  nom- 
inally operated  about  the  0.4  amps  and  2.0  volts  operating  point,  from  [Ref  12]  their 
operation  curves  and  technical  data  is  available. 

3.  Heaters 

As  a  means  of  degassing  the  FC-72  liquid  and  for  bulk  temperature  control,  two 
strip  heaters  of  125mm  x  1 1  mm  were  employed.  These  strip  heaters  were  secured  to  the 
bottom  interior  surface  of  the  chamber  with  epoxy,  and  were  orientated  lengthwise  in  the 
chamber  such  that  one  heater  strip  will  be  in  front  of  and  one  heater  strip  behind  the 
insert  board.  To  raise  the  bulk  fluid  temperature  to  saturation  temperature  and  for 
degassing  purposed,  both  heaters  are  on  line  and  operated  at  1.8  amps  and  20  volts. 
Once  degassing  is  complete,  the  forward  strip  heater  is  electrically  isolated  and  the  rear 
heater  is  run  at  0.9  amps  and  20  volts.  This  will  maintain  the  forward  portion  of  ex- 
perimental chamber  at  near  saturation  temperature  with  minimum  stratification  of  fluid, 
and  maintain  liquid  in  a  degassed  state  since  the  fluid  is  still  boiling  from  the  heater 
strip's  surface. 

4.  Thermocouples 

A  total  of  Six  copper-constantan  type  thermocouples  with  a  wire  diameter  of 
0.05  inches  (0.127  mm)  were  utilized  for  temperature  measurements.  Two  thermocouples 
were  bonded  to  the  underside  of  the  aluminum  condensor  plate  utilizing  OMEGA  Bond 
type  101  thermally  conductive  epoxy,  four  thermocouples  were  installed  on  the  insert 
board,  and  one  was  used  for  measuring  the  ambient  air  temperature. 

5.  Insert  Board 

An  insert  board  was  designed  to  mount  into  two  0.250  inch  wide  grooves  that 
were  milled  into  the  chamber  side  walls.     Onto  this  insert  board  are  mounted  four 


thermocouples  to  measure  bulk  fluid  temperature,  five  platinum  wires  of  nominally  0.002 
inch  (0.05mm)  diameter  by  4.0  inch  long  for  conducting  the  actual  boiling  experiments. 
Five  springs  are  used  to  maintain  the  platinum  wire  tension  so  no  sagging  develops  due 
to  thermal  expansion  of  the  wires.  In  Appendix  A  the  platinum  wire  calibrations  and 
specifications  are  shown  in  detail. 

B.     INSTRUMENTATION 

In  the  design  and  assembly  of  this  experimental  apparatus,  five  channels  were  wired 
and  instrumentation  provided  to  measure  voltage  drops  at  select  points.  Each  of  these 
channels  had  a  2.0  ohm  precision  resistor  wired  in  series  with  the  platinum  wire,  where 
the  voltage  drop  across  the  resistor  was  used  to  obtain  the  current  flow  through  the 
platinum  wire.  The  voltage  drop  measured  across  the  platinum  wire  is  used  to  obtain 
the  wire  resistance,  power  dissipated  and  the  heat  flux  of  the  platinum  wire.  The 
platinum  wire  temperature  is  obtained  from  its  resistance  vs  temperature  calibration 
curve  obtained  prior  to  experimental  runs.  Instrumentation  was  additionally  used  to 
obtain  the  bulk  fluid  temperature,  which  will  be  required  for  calculating  fluid  properties. 

1.     Platinum  Wire 

Each  of  the  five  platinum  wires  were  run  in  series  with  a  2  ohm  (+  1%)  precision 
resistor  of  25  watt  power  rating.  The  following  table  gives  the  platinum  wire  and  resistor 
combinations  which  were  used. 

Table   I.     MEASURED  RESISTANCE  VALUES  FOR  2  OHM  RESISTORS 


Wire  Chan- 
nel 

Measured  Resistance,  20  Sample 
average 

20  Sample  Std.  Dev. 

1 

I.9994Q 

0.0009Q 

2 

2.0076Q 

0.000512 

3 

1.9998Q 

0.0003Q 

4 

2.0018Q 

0.001 3Q 

5 

2.0025Q 

0.0008Q 

A  2.0  ohm  precision  resistor  was  chosen  over  the  previous  design's  0.1  ohms 
precision  resistor  in  order  to  reduce  the  calculated  currents  uncertainty,  and  there  by 
reduce  the  calculated  platinum  wire  temperature  uncertainty,  [Ref.  11:  p.  25]. 


2.  Power  Supplies 

The  following  power  supplies  were  utilized  and  their  primary  function  listed. 

1.  Hewlett  Packard  6214C  0-lOvolts/O-l  amps:    Used  to  power   the  thermal  electric 
coolers. 

2.  Hewlett  Packard  6214C  l-10volts  /  0-1  amps:  Used  to  power  the  constant  heat  flux 
platinum  wire. 

3.  WP711  0-40  volts  /  0-1  amps  :  Used  to  power  the  varying  heat  flux  platinum  wire. 

4.  Hewlett  Packard  6286A  0-20volts  /  0-10  amps  :  Used  to  power  Heater  Strips  for 
degassing  and  bulk  temperature  maintenance. 

5.  Hewlett  Packard  6289A  0-40  volts  /  0-1.5  amps  :  Used  to  power  an  auxiliary  illu- 
mination system  or  an  augmenting  cooling  fan  for  the  thermal  electric  coolers. 

3.  Acquisition  Unit 

To  obtain  the  necessary  measurements,  a  Hewlett  Packard  3852A  Data  Acqui- 
sition System  was  employed,  utilizing  a  Hewlett  Packard  300  series  desktop  computer. 
The  HP3852A  utilized  the  following  utility  modules  inserted  into  it's  chassis. 

1.  Hewlett  Packard  Model  44701A  Digital  Integrating  Voltmeter 

2.  Hewlett  Packard  Model  44705A  20  channel  relay  multiplexer  for  measuring  re- 
quired voltage  drops. 

3.  Hewlett  Packard  Model  447 13A  24  channel  relay  multiplexer  with  electronic  cold 
junction  for  direct  conversion  of  thermal  couple  emf  to  temperature. 

The  HP300  desktop  computer  was  used  to  control  the  HP3852A  Data  acquisition  sys- 
tem. Additionally,  through  specifically  written  programs,  the  computer  can  accomplish 
the  required  data  reductions,  (See  Appendix  D). 
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Figure  2.      Overall  experimental  settup 


Figure  3.      Experimental  Chamber  in  Detail 
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Figure  4.     Side  View  of  Experimental  Apparatus 
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Figure  5.      Rear  View  of  Experimental  Chamber 
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Figure  6.      Insert  Board 
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Figure  7.      schematic  Drawing  of  Insert  Board 


14 


Figure  8.      Isometric  View  of  Apparatus 
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Figure  9.      Condensor  Surface  Details 
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III.     EXPERIMENTAL  PROCEDURE 

A.     PREPARATIONS 

Careful  preparations  were  required  to  ensure  that  the  data  obtained  was  of  high 
quality.  The  most  notable  requirement  being  the  necessity  of  thoughly  degassing  the 
FC-72  liquid,  [Ref,  3:  p.  3]. 

1.  Normal  Warm-Up  Procedure 

In  carrying  out  the  experimental  runs  the  following  procedures  were  followed. 

1.  Power-up  the  HP3852A  Data  Acquisition  Unit  to  start  the  warm-up  period  of  one 
hour  minimum.  [Ref.  13:  p.  586] 

2.  Ensure  the  liquid  level  is  above  Wire  5,  but  below  the  condensate  drip  tray. 

3.  Set  power  supply  of  thermoelectric  coolers  to  1.3  amps  and  2.0  volts. 

4.  Ensure  the  experimental  chamber's  insulation  jacket  is  in  place  and  set  power  sup- 
ply to  strip  heaters  to  approximately  1.8  amps  and  20  volts. 

5.  Verify  the  0.25  inch  diameter  vent  line  is  unobstructed. 

6.  Power  on  the  HP300  series  desktop  computer,  and  load  data  acquisition  program 
as  contained  in  Appendix  D. 

7.  Once  fluid  has  reached  saturation  and  is  boiling  vigourisly,  hold  for  a  minimum  of 
one  hour,  typically  two  hours  was  achieved,  to  permit  adequate  degassing  of  fluid. 

In  actual  practice,  the  time  required  to  raise  the  bulk  temperature  of  the  fluid  to  satu- 
ration temperature  of  56  °C  was  nominally  two  hours. 

2.  Post  Warm-up  Procedure 

Once  the  required  degassing  has  been  accomplished,  power  is  temporarily  se- 
cured to  the  strip  heaters.  Once  power  is  secured  the  forward  strip  heater  is  electrically 
isolated.  The  rear  strip  heater  is  powered  at  0.9  amps  and  20  volts,  and  the  thermo 
electric  coolers  are  secured. 

The  power  setting  of  the  strip  heaters  permits  vigorous  boiling  from  it's  surface. 
This  serves  to  keep  the  fluid  degassed  and  the  bulk  fluid  temperature  near  saturation  and 
near  isothermal.  By  keeping  the  fluid  level  below  the  insert  board's  upper  edge,  the  fluid 
forward  of  the  insert  board  remains  undisturbed  by  the  boiling  action  occurring  in  the 
rear  chamber. 
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B.  DATA  ACQUISITION  PROCEDURE 

Upon  completion  of  the  previous  preparatory  procedures,  the  experimental  runs  are 
ready  to  begin.  The  procedure  followed  for  an  experimental  run  is  as  follows: 

1 .  Plug  power  supply  into  power  receptacles  for  wire  1 .  Power  on  and  ensure  current 
limiter  is  set  to  max  value  of  1.0  amps. 

2.  Plug  HP6214C  power  supply  into  power  receptacles  of  desired  constant  heat  flux 
wire. 

3.  Set  power  supply  of  step  2  to  proper  output  to  give  desired  heat  flux  from  the  wire. 

4.  Start  running  the  Data  Acquisition  Program,  which  is  programmed  to  measure  the 
following  data  from  the  desired  wires  at  two  minute  intervals. 

•  Precision  Resistor  Voltage  drop,  up  to  five  channels. 

•  Platinum  Wire  Voltage  drop,  up  to  five  channels. 

•  Four  bulk  fluid  temperature  measurements 

•  Two  condenser  surface  temperature  measurements. 

•  One  ambient  air  temperature  measurement. 

5.  After  each  data  interval,  increase  the  voltage  to  the  varying  heat  flux  wire  by  0.2 
Volt  increments.  This  was  continued  until  the  on  line  power  supply  was  at  2.3  volts 
of  outputs,  the  increments  were  then  reduced  to  0.1  volts  until  nucleation  occurred. 
The  increments  were  then  returned  to  0.2  volts  until  a  heat  flux  of  approximately 
200,000  watts/m2  were  achieved.  At  any  higher  heat  flux  the  wire  will  go  into  the 
film  boiling  regime. 

6.  Reduce  the  power  to  the  varying  heat  flux  wire  at  0.3  Volts  increment  to  bring  the 
wire  back  to  a  no  power  condition. 

7.  Repeat  steps  3  thru  6  for  all  other  wire  and  heat  flux  configurations. 

C.  DATA  REDUCTION 

Data  that  was  obtained  from  the  actual  experimental  apparatus  were  two  voltage 
measurements  per  active  platinum  wire  channel,  and  six  temperature  measurements. 
The  process  used  for  data  reduction  is  a  combination  of  computer  programs  and 
graphical  outputs  for  numerical  presentation  and  visual  presentations  respectively. 

1.     Numerical  Data 

The  numerical  data  values  required  to  be  processed  are: 

1.  Fluid  Bulk  Temperature 

T^-tt+V^ra  0..) 

and  where,  T,  thru  T4  are  the  thermocouple  measurements  of  the  bulk  fluid. 

2.  Current  flow  through  the  Platinum  Wire,  where 
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hire  =  T^  (Amps) 

3.   Platinum  Wire's  electrical  resistance  value,  where 

V, 


R 


PtWire 


PrlVire 

*  ur'tro 


(Ohms) 


4.  Platinum  Wire  surface  temperature,  where 


TjUr  —  ^PtWire  +  T0(°C) 


(3.2) 


(3.3) 


(3.4) 


where, 

T„r=  the  actual  wires  surface  temperature 

£  is  the  slope  of  the  calibration  curve 

T„  is  the  calibration  curves  y-axis  intercept 


5.  Fluid  Film  Temperature,  where 


T«/m  — 


Trur  +  ^butk 


Jllm  2.0 

6.  Dissipated  Heat  Flux  of  platinum  wire,  where 


q  = 


IvWre  x  V PtWire    ,   Watts 


Area 


(3.5) 


(3.6) 


Surf 


m 


where,  Areaw  =  n  x  lengthw/r,  x  D, 


The  above  numerical  data  were  all  derived  from  the  raw  data,  and  previously 
determined  calibration  curves.  The  following  table  provides  physical  dimensions  and 
temperature  calibration  coefficients  of  the  platinum  wires. 


Table  2.     PLATINUM   WIRE   PHYSICAL   PARAMETERS    FOR   EQUATION 
(3.4)  AND  EQUATION  (3.5) 

Wire 

e 

T0 

length 

diameter 

°C/Q 

°C 

mm 

mm 

1 

53.902 

-254.470 

99.22 

0.05 

2 

53.666 

-253.009 

99.22 

0.05 

3 

53.141 

-252.729 

100.01 

0.05 

4 

53.484 

-252.848 

98.43 

0.05 

5 

53.086 

-260.049 

100.01 

0.05 

To  see  in  detail  the  methodology  used  to  obtain  the  numerical  values  presented 
in  the  preceding  table,  proceed  to  Appendix  A  'Platinum  Wire  Calibration'. 
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2.     Fluid  Property  Data 

The  Fluid  Properties  to  be  determined  are  as  follows: 

1.  Thermal  Conductivity  of  fluid,  as  obtained  from  Figure  5  of  the  3M  Corporation 
Flourinert  Product  Manual  [Ref.  3:  p.    15]  have  been  determined  to  be, 

0.6033  -0.001 15  *T>,        watts, 
k  = K (-^-)  (3.7) 

2.  Thermal  Expansion  Coefficient,  from  table  4B  in  the  Product  Manual  [Ref.  3:  p.  10] 
is  given  as 

n  0-00261  ,    1    v  n„v 

p    '  1.740-0.00251  xTy^  l  °C  ]  K      ' 

3.  Kinematic  Viscosity  of  the  Fluid  ,  as  obtained  from  Figure  (3)  in  the  Product  Man- 
ual [Ref  3:  p.  13]  and  determining  an  exponential  curve  fit  yields 

v=1.203952E-8  x  exp(l,058.4109°- L.    _0     )(-^-)         (3.9) 

lbuik  +  213-15  c        5 

4.  Specific  Heat  of  the  Fluid  for  all  the  Flourinert  Fluids  are  obtained  from  Figure  (4) 
of  the  Product  Manual  [Ref  3  :  p.  14]  is, 

Cp  =  (0.24111  +  3.70374E-  4  x  Tbulk)  x  4186(t4tt)  (3.10) 

Kg    v^ 

5.  Liquid  Density  of  FC-72,  from  expression  in  table  4B  and  constants  from  table  4C 
in  the  Product  Manual  [Ref  3:  p  10]  is  given  by 

ke 
p  =  (1.740  -0.00261  xTM)x  1000(— ^)  (3.11) 

m 


6.  Thermal  Diffusivity  of  the  liquid  is  given  by, 

k 

pxCp 


«  =  — ^F~  (3-12) 


7.  Prandlt  Number  of  the  Fluid  is  given  by 


Pr  =  x  (3.13) 

8.  Ralyeigh  Number  of  Fluid  based  on  wire  diameter  is  given  by, 

gxflx  (T„.r  —  Ihuw)  x  D 
Ra  =  v  x  «        (3-14) 

9.  Correlation  for  the  Nusselt  Number  was  obtained  by  using  the  correlation  of  Kuehn 
and  Goldstein  [Ref  10:  p.  1128]  given  by 

Nu^  = 2  2 

m*  [(0.518xRai/4x[l+(-0^)3/5^5/12)15  +  (0.1xRal/3)15]1/,5  ] 
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10.  Theoretical  Heat  Flux  as  obtained  from  using  the  ND  value  obtained  from  equation 
(3.15). 

k  x  Nug  x  (Tsur  -  Tbu[k)     watts  >  „  ,,, 

u  m 

With  all  the  required  numerical  data  and  fluid  properties  obtained,  it  is  now 
necessary  to  present  the  numerical  information  in  a  graphical  representation  to  see  the 
relevancy  of  the  results  obtained.  This  will  become  apparent  with  the  presentation  of 
results  in  Chapter  IV  'Results  and  Discussion'. 


21 


IV.     RESULTS  AND  DISCUSSION 

The  data  obtained  can  be  grouped  for  discussion  as  follows: 

1.  Individually  powered  wires 

Saturated  Conditions 
Sub-Cooled  Conditions 

2.  Bubble  Pumping  Effect  Runs 

Saturated  Conditions 
Sub-Cooled  Conditions 

A.     INDIVIDUALLY  POWERED  WIRES 

Before  analyzing  the  pumping  effect  of  the  bubbles  on  the  upper  wires,  the  individ- 
ual wire  boiling  experiments  will  be  investigated.  All  of  the  individual  powered  wire  runs 
shown  here  were  made  at  or  near  saturated  fluid  conditions  with  two  exceptions.  The 
following  table  will  provide  a  listing  of  single  wire  runs. 


Table  3.     INDIVIDUAL  WIRE  RUNS 


Figure 


Wire 


*  bulk 


Figure  10 


1 


55.5°  C 


Figure  11 


1  -  several  runs  overlayed 


At  or  near  saturation 


Figure  12 


1 


38°  C 


Figure  13 


23.0°  C 


Figure  14 


56.0°  C 


Figure  15 


56.0°  C 


Figure  16 


56.0°  C 


Figure  17 


all  the  above 


see  individual  wires 


In  all  cases  conditions  were  under  atmospheric  pressure  with  0.05mm  diameter  platinum 
wires. 

1.    Natural  Convection 

Figure  10  is  a  typical  boiling  curve  for  FC-72  with  the  Kuehn  and  Goldstein 
natural  convection  region  shown  along  with  the  experimental  points.  The  nucleate 
boiling  zone  starts  at  boiling  incipience,  which  is  evident  from  the  sudden  drop  of  wall 
temperature,  a  characteristic  of  a  highly  wetting  fluid.  Figure  His  a  compilation  of  wire 
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1  data,  at  or  near  saturation,  obtained  over  a  period  of  time.  As  can  be  seen,  with  the 
exception  of  the  boiling  incipience  point,  the  data  is  virtually  super  imposable  over  each 
other.  This  reproducability  of  data  over  time  indicates  no  significant  aging  has  occured 
on  wire  1,  and  more  importantly  no  change  in  calibration. 

In  fact  in  the  natural  convection  region  Figure  10  thru  Figure  16  are  all  essen- 
tially superimposable  over  each  other.  Figure  17  is  an  overlay  of  all  the  previous  fig- 
ures. Therefore,  it  appears  that  whether  the  fluid  is  Saturated  or  Sub-Cooled,  the  fluid 
bulk  temperature  appears  to  have  little  effect  on  the  natural  convection  heat  transfer 
rate. 

There  was  generally  very  good  agreement  between  the  natural  convection  data 
and  the  Kuehn  and  Goldstein  correlation  [Ref  10]  (equation  3.15).  The  only  deviation 
existing  between  the  figures  being  the  total  temperature  excursion  to  incipience,  which 
ranged  from  27°  C  to  38°  C. 

2.  Nucleate  Boiling 

As  previously  stated,  the  nucleate  boiling  zone  starts  at  boiling  incipience.  Us- 
ing wire  1  as  a  standard  and  Figure  lOas  base  line  boiling  characteristic  data  a  compar- 
ison to  other  wires  at  saturation  condition,  and  a  comparison  to  subcooling  effects  on 
wire's  1  boiling  can  be  made.  As  is  evident  from  Figure  17,  the  boiling  characteristic 
of  all  the  wires  at  saturation  condition  tend  to  overlap  each  other.  The  differences  aris- 
ing tending  to  be  caused  by  individual  wire  placement  in  the  fluid  .  However,  the  effect 
of  subcooling  does  result  in  a  decrease  in  the  boiling  curve  overshoot  at  point  of 
incipience.  This  decrease  in  the  boiling  curve  overshoot  results  in  a  net  decrease  in  the 
superheat  required  to  initiate  nucleate  boiling. 

3.  Departure  from  Nucleate  Boiling 

The  Critical  Heat  Flux  (CHF)  for  FC-72  at  saturated  condition  was  calculated 
by  Danielson  et  all  [Ref  9]  to  be  approximately  203,000  Watts/m2.  Experimental  value 
for  CHF  ranged  from  198,000  Watts/m2  to  220,000  Watts/m2,  a  spread  of  -2.5%  to 
8.4%.  A  CHF  for  a  subcooled  condition  of  T4l//A=38°  C  was  found  experimentally  to 
be  approximately  400,000  watts/m2.  A  CHF  for  Tbulk=  23  °  C  was  not  obtained  during 
the  experimental  runs. 
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Figure  10.      Wirel  boiling  curve  Tbulk  =  55.5°C 
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Figure  11.      Wire  1  Boiling  Curves  over  time 


25 


100000 


CN 

2 


§ 


X 

3 


< 

UJ 

X 


10000 


1000- 


^ 


III 

ft 


../.. 


./. 


WIRE  1  BOILING  CURVE 


-O-O- 

a 


o   ;a 


o      ° 

_a 


.a. j.. a. 

;a— . 


"7a 


-9   5 


-o  4 


-•    3 


-«    2 

-°  CD 


LEGEND 
a  =  HEAT  FLUX1  INCREASE 
o  =  HEAT  FLUX1  DECREASE 
GOLDSTEIN  CORRELATION 


10 


20 


—r- 
25 


- 1 — 
30 


35 


TSURF  -  TBULK  (C) 


40 


Figure  12.      Wire  1  boiling  curve  T4u/A  =  38°C 
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Figure  13.      Wire  1  Boiling  Curve  Tbulk  =  23°C 
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Figure  14.      Wire  2  boiling  curve  Tbulk  =  56°C 
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Figure  15.      Wire  3  Boiling  Curve  7bulk  =  56°C 
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Figure  16.      Wire  4  boiling  curve  Tbulk  =  56°C 
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Figure  17.      Overlay  of  previous  curves 
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B.      BUBBLE  PUMPING  EFFECT  ON  HEAT  TRANSFER  RATE 

In  the  upcoming  discussion  a  brief  explanation  of  the  figures  are  required.  In  each 
case  for  these  experiments,  heater  wire  number  1  was  used  to  generate  the  bubbles  for 
examining  the  pumping  effects  upon  the  upper  wires.  The  figure  insert  is  used  to  show 
which  wires  are  being  utilized.  Only  select  data  points  were  plotted  for  the  passive 
(constant  heat  flux  wire)  wire,  and  those  points  correspond  in  time  to  the  identically 
numbered  data  point  on  wire  l's  boiling  curve.  The  data  points  were  chosen  since  they 
defme  regions  of  wire  l's  boiling  curve  as  shown  in  the  following  table. 

Table  4.     DATA  POINT  DESCRIPTION 


5  point  figure 

6  point  figure 

Region  Description 

1 

1 

Natural  Convection  Start 

2 

2 

Point  of  Incipience 

3 

Partial  Boiling  on  wire 

3 

4 

Full  nucleate  Boiling 

4 

5 

Max  Nucleate  Boiling 

5 

6 

End  of  run 

1.     Passive  Wire  at  Low  Heat  Flux 

In  these  experimental  runs  the  heat  flux  on  the  passive  wire  was  set  to  approx- 
imately 8500  watts/m2.  At  this  low  heat  flux  it  was  expected  that  little  to  no  active 
nucleation  will  occur  on  the  passive  wire,  even  after  wire  1  is  in  nucleate  boiling. 
Therefore,  all  the  enhancement  of  the  upper  wire  will  be  due  to  an  addition  of  a  force 
convection  term  to  the  heat  transfer  rate. 

In  comparing  Figure  18  to  Figure  19the  following  can  be  inferred: 

1.  Total  reduction  of  wall  temperature  on  the  passive  wires  are  comparable.    Wire  2 
having  a  3.7°C  and  Wire  4  having  a  3.0°  C  drop  in  wall  temperature. 

2.  A  definite  effective  distance  exist  over  which  wire's  1  influence  is  imparted  upon  the 
upper  wires,  for  example: 

Wire  2  and  wire  4  showed  a  very  small  change  in  wall  temperature  at  point  2  (point 
of  incipience). 

After  incipient  boiling  has  occurred  on  wire  1,  wire  2  showed  a  greater  drop  in  wall 
surface  temperature  then  wire  4.  Approx.  2.5°  C  drop  for  wire  2  and  0.1°  C  drop 
for  wire  4. 
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3.  At  point  4,  the  point  of  maximum  temperature  in  nucleate  boiling,  both  wire  2  and 
wire  4  showed  a  comparable  surface  temperature,  approximately  2°  C  above  satu- 
ration. 

As  can  be  seen  from  the  figures,  some  nucleation  sites  have  been  activated  on  wire  2  and 
wire  4,  with  point  5  showing  a  lower  wall  temperature  then  point  1. 

2.  Passive  Wire  at  Medium  Heat  Flux 

In  these  experimental  runs  the  heat  flux  on  the  passive  wire  was  set  to  approx- 
imately 20,500  watts/m2.  At  this  heat  flux  it  was  expected  that  active  nucleation  will 
occur  on  the  passive  wire  after  wire  1  is  nucleate  in  boiling. 

In  comparing  Figure  20  to  Figure  21 
the  following  can  be  inferred: 

1.  The  total  reduction  of  wall  temperature  on  the  passive  wire  ,  points  1  to  points  4, 
are  comparable  at  7.1°  C  for  wire  2,  and  6.3°  C  for  wire  4. 

2.  The  spatial  effect  is  evident  on  the  influence  of  wire  1  on  the  upper  wires,  for  ex- 
ample: 

Between  point  1  and  point  2,  wire  2  shows  a  greater  reduction  in  wall  surface  tem- 
perature than  is  evident  in  wire  4. 

At  point  3,  wire  2  shows  a  much  greater  drop  in  surface  temperature  than  wire  4: 
4°  C  compared  to  2°  C.  The  difference  being  the  amount  of 

3.  At  point  4,  maximum  nucleate  boiling  from  wire  1,  both  wire  2  and  wire  4  show 
comparable  values  for  surface  temperature,  approximately 

1.3°  C. 

As  can  be  seen  from  the  figures,  wire  1  has  activated  previously  inactive  nucleation  sites 
on  wire  2  and  wire  4.  This  is  evident  from  point  5  being  at  a  lower  surface  temperature 
than  point  3. 

3.  Passive  Wire  at  High  Heat  Flux 

In  these  experimental  runs  the  heat  flux  on  the  passive  wire  was  set  to  40,000 
Watts/m2  on  wire2  and  35,000  Watts/m2  on  Wire  4.  At  this  heat  flux  it  was  expected 
that  a  large  number  of  inactive  nucleation  sites  will  be  activated  after  wire  1  is  nucleate 
in  boiling. 

In  comparing  Figure  22  and  Figure  23  it  is  evident  that  the  information  ob- 
tainable from  the  figures  is  comparable  to  the  information  for  the  medium  heat  flux  case. 
However,  the  wall  temperature  reductions  on  the  passive  wires  are  much  greater  as 
compared  to  the  Medium  Heat  Flux  case,  in  particular  the  following: 

1.  In  comparing  the  surface  temperature  reduction  from  point  1  to  point  2,  wire  2  has 
dropped  2.5°  C  while  wire  4  had  almost  a  zero  drop  of  temperature.    This  appears 
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to  be  due  to  the  relative  locations  of  the  passive  wires  with  respect  to  wire  1,  and 
is  especially  apparent  at  the  higher  heat  fluxes. 

2.  Many  inactive  nucleation  sites  have  been  activated  by  wire  l's  nucleation  bubbles. 
As  can  be  seen  by  comparing  point  1  to  point  5(6),  a  difference  on  Wire  2  of  13° 
C  and  on  Wire  4  of  7°  C  is  apparent. 

4.  Wire  1  Passive  with  Wire  2  varying 

In  these  experimental  runs,  wire  1  was  set  to  150,000  Watts/m2  so  it  was  boiling 
vigorously,  wire  2  was  actively  varied  over  a  large  heat  flux  range  as  can  be  seen  from 
Figure  24.  On  Figure  24  a  previously  obtained  boiling  curve  of  Wire  2  was  overlayed 
onto  it.  Upon  comparing  the  two  curves  the  most  striking  feature  is  the  disappearance 
of  the  Hysterises  Loop  as  is  evident  in  the  overlayed  boiling  curve  data.  In  this  ar- 
rangement, the  boiling  curve  overshoot  has  been  removed  as  a  design  problem,  bubbles 
from  Wire  1  has  removed  a  major  design  consideration. 

5.  Subcooled 

In  these  experiments,  the  bulk  fluid  was  degassed  and  then  subcooled  to  a  tem- 
perature of  38°  C.  In  these  experimental  runs,  wire  2  was  set  at  35,000  watts/m2  and 
53,000  watts/m2,  as  shown  in  Figure  25  and  Figure  26  respectively.  In  subcooling,  no 
significant  effects  were  notable  beyond  those  seen  in  the  saturated  runs.  The  only  effect 
that  was  apparent  was  the  ability  to  run  wire  2  at  a  greater  heat  flux  with  out  incipent 
boiling  occuring  on  the  wire. 
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Figure  18.      Wire2  at  8500  W/  m2  Tbulk  =  55.5°C 
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Figure  19.      Wire4  at  8500  W/  m2  Tbulk  -  55.6° C 
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Figure  20.      Wire  2  at  20,500  W/  m2  lbulk  =  55.9degreeC 
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Figure  21.      Wire  4  at  20,500  W/  m2  T4urt  =  55.6°C 
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Figure  22.      Wire  2  at  40,000  W/  mJ  Tbutk  =  54. 6°  C 
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WIRE  4  AT  35,000  WATT/M2 

NEAR  SATURATION  Tbulk  =  55.0°Celisus 


100000 


en 

5 


x 

3 


10000 


1000 


TSURF  -  TSAT  (C) 


Figure  23.      Wire  4  at  35,000  W/  m2  Tbulk  =  55°C 
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Figure  24.      Wire  1  passive,  Wire  2  active  Tbulk  =  56°  C 
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Figure  25.      Wire  2  at  35,000  W/  m2  Tbulk  =  38°C 
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WIRE  2  AT  53,000  WATT/M2 
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Figure  26.      Wire  2  at  53,000  W/  m2  Tbulk  =  38°C 
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V.     CONCLUSIONS 

The  study  conducted  in  this  thesis  was  to  investigate  the  effects  of  the  wake  from  a 
bubble  generating  wire  on  wires  located  downstream  within  the  wake.  Experimental 
runs  were  conducted  at  saturated  condition  and  subcooled  conditions  in  which  the  fol- 
lowing were  found: 

A.  SINGLE  WIRE  BOILING 

1.  The  effects  of  subcooling  on  the  natural  convective  heat  transfer  from  the  wire  was 
negligible. 

2.  The  effects  of  subcooling  tended  to  increase  the  heat  flux  required  for  incipient 
boiling  to  initiate. 

3.  With  subcooling,  boiling  hystersis  is  still  observed,  but  the  hystersis  magnitude  de- 
creases as  subcooling  increases. 

4.  Subcooling  results  in  a  reduction  in  wall  temperature  after  the  start  of  nucleate 
boiling  as  compared  to  saturated  cases. 

B.  PASSIVE  WIRE  ANALYSIS 

Heat  Transfer  enhancement  effects  from  Wire  1  upon  the  downstream  wires  has 
many  factors  within  it.  The  following  factors  have  been  identified: 

1.  An  enhancement  to  convection  resulting  from  buoyancy  induced  fluid  flow  off  of 
wire  1  was  seen  to  exist.  This  effect  is  very  spatially  dependent  with  the  effects 
most  pronounced  on  wire  2  and  negligible  on  wire  4.  This  corresponds  to  point  1 
and  point  2  on  the  figures. 

2.  The  sudden  drop  of  wall  temperature  between  point  2  and  point  3/4  was  seen  to 
be  the  results  of  the  following: 

Nucleation  bubbles  from  wire  1  activating  the  inactive  nucleation  sites  on  the  upper 
wires  resulting  in  boiling  heat  transfer  off  of  the  upper  wires. 

An  additional  force  convection  enhancement  due  to  the  nucleation  bubbles  from 
wire  1  breaking  up  the  thermal  boundary  layer  about  the  upper  wires  is  believed  to 
exist. 

3.  The  drop  in  temperature  on  the  upper  wires  corresponding  with  points  3/4  to  point 
4/5  appears  to  be  due  purely  to  an  increase  in  force  convection  due  to  the  high 
density  of  bubbles  off  of  wire  1. 

Additional  conclusions  reached  are. 

1.  A  constant  source  of  bubble  generation  from  below  is  an  effective  way  to  eliminate 
the  Boiling  Curve  Overshoot  and  Hysteresis  loop  associated  with  dielectric  fluids. 
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At  low  heat  fluxes,  proximity  to  Wire  1  is  important  for  heat  transfer  enhancement 
to  occur.   This  dependency  reduces  as  heat  flux  increase  on  the  passive  wire. 
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APPENDIX  A.     CALIBRATION 

A.     PLATINUM  WIRE 

1.    Surface  Micro-Structure 

Prior  to  the  experiments,  a  length  of  0.05mm  platinum  wire  was  analyzed  using 
a  Scanning  Electron  Microscope  (SEM).  The  length  of  wire  analyzed  was  not  used  in 
any  experiments,  but  it  did  come  from  the  same  spool  of  wire  as  the  experimental  lengths 
of  wire,  therefore,  it  should  be  a  representative  sample.  In  analyzing  this  sample  the 
SEM  brought  out  several  important  features: 

1.  From  Figure  27  a  wire  diameter  of  0.052mm  can  be  determined.  Also,  the  surface 
is  relatively  unmarked  at  this  magnification  of  1490x;  little  scoring  on  the  wire  and 
a  low  density  of  surface  pits  are  evident. 

2.  In  Figure  28,  an  increase  in  magnification  to  4740x  within  the  field  of  Figure  27, 
a  more  detailed  view  of  the  larger  surface  defects  are  evident.  Maximum  defect  size 
is  approx  2^m  in  length  and  O.Sfim  in  width.  A  surface  defect  of  that  size  would 
be  too  large  to  be  a  possible  nucleation  site  ,  according  to  the  calculations  done  in 
Ref.  4of0.1/tm. 


3. 


In  Figure  29,  a  further  increase  in  magnification  to  13,200x  within  the  field  of 
Figure  27,  pits  of  approximately  0.8jmi  max  down  to  0.3^m  are  seen.  It  will  be 
these  surface  defects  that  boiling  will  initiate  from.  The  density  of  these  pits  along 
the  length  of  the  wire  analyzed  was  not  great,  this  implies  that  the  incipient  super- 
heat will  be  large. 


Figure  27.      Surface  micro  geometry:     Sample  Platinum  Wire  (1490x) 
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Figure  28.     Surface  micro  geometry:    Sample  Platinum  Wire  (4740x) 


Figure  29.      Surface  micro  geometry:     Sample  Platinum  Wire  (13,200x) 
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2.     Desirability  for  Using  Platinum 

Platinum  was  chosen  for  several  reasons: 

1.  It  has  a  highly  linear  response  of  Temperature  vs  Resistance. 

2.  It  has  a  large  AR  for  a  given  AT. 

In  each  aspect,  platinum  performed  very  well.  As  can  be  seen  in  the  calibration  curves 
of  Figure  31  thru  Figure  40  response  is  linear,  with  a  0.1Q.  change  in  resistance  for  ev- 
ery 5°  C  change  in  temperature. 

B.  BACKGROUND 

An  important  step  in  minimizing  experimental  data  uncertainty  is  an  accurate  cali- 
bration of  the  platinum  wires  used.  In  obtaining  the  calibration  curves  for  the  wires, 
three  calibration  runs  were  performed. 

1.  An  initial  calibration  in  a  heated  calibration  bath  with  a  platinum  resistance 
thermometer. 

2.  An  in  place  calibration  of  wires  in  the  test  chamber,  with  the  four  installed 
thermocouples  used  for  temperature  measurements,  was  conducted  to  ensure  the 
calibration  curves  from  the  calibration  bath  were  not  altered  from  handling  the  wire 
insert  board. 

3.  A  final  in  place  calibration  of  wires  in  the  test  chamber, 

with  the  four  installed  thermocouples  used  for  temperature  measurements,  was 
conducted  to  ensure  the  calibration  curves  used  for  data  acquisition  were  not  al- 
tered due  to  any  reason. 

The  second  calibration  run  was  required  to  ensure  no  distortion  occurred  to  the  wires 
when  the  insert  board,  on  which  the  wires  are  mounted,  was  removed  from  the  calibrated 
bath  and  installed  into  the  experimental  test  chamber. 

C.  INITIAL  CALIBRATION 

The  initial  calibration  was  conducted  using  the  following  equipment. 

1.  Rosemount  Engineering  Co.  Model  913A  calibration  bath  with  Etheylene  Glycol 
as  the  working  fluid. 

2.  Rosemount  Engineering  Co.  Model  923B  power  supply,  controlling  calibration 
bath  temperature. 

3.  Rosemount  Engineering  Co.  Model  920A  commutating  bridge,  used  to  measure  an 
resistance  to  0.001Q  accuracy. 

4.  Rosemount  Engineering  Co.  Model  162C  S/N  985  Platinum  Temperature  Probe. 
Used  to  measure  calibration  bath  temperature  through  the  commutating  bridge. 

5.  Hewlett  Packard  Data  Acquisition  System.  Used  to  measure  platinum  wire  resist- 
ance using  a  4-wire  ohm  measurement  technique. 
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1.  Calibration  Heating  Bath 

The  present  investigation  was  conducted  using  an  externally  heated  calibration 
bath  with  etheylene  glycol  as  the  working  fluid.  This  fluid  was  circulated  by  a  centrifugal 
pump  to  ensure  a  uniform  bulk  temperature  throughout  the  bath.  Item  (1)  and  item  (2) 
are  utilized  . 

2.  Temperature  Measurements 

Temperature  measurements  for  calibration  were  conducted  using  a  high  sensi- 
tivity commutating  bridge  to  measure  the  resistance  of  a  platinum  resistance  temper- 
ature probe.  This  platinum  temperature  probe  had  been  previously  calibrated  to  a 
standard  that  is  traceable  to  the  National  Bureau  of  Standards. 

3.  Resistance  Measurements 

In  calibrating  the  platinum  wire  the  resistance  of  the  wire  was  measured  with  the 
same  Data  Acquisition  Unit  Components  as  would  be  used  in  the  experimental  runs. 
This  was  to  ensure  the  uncertainty  of  measured  resistance  would  be  the  same  as  between 
the  calibration  run  and  experimental  runs. 

The  Data  Acquisition  Unit  Component  used  are: 

1.  HP3852A  Control  Unit 

2.  HP44701A  Integrating  Voltmeter  card 

3.  HP44705A  20  channel  Multiplexer  card 

4.  HP300  series  desk  top  computer 

5.  HP44713A  24  channel  multiplexer  with  electronic  cold  junction  for  thermo  couple 
measurements. 

4.  Calibration  Procedure 

By  using  the  previously  described  components  for  measuring  temperature  and 
resistance,  the  following  procedure  was  used  to  obtain  the  system's  calibration  data: 

1.  Set  calibration  bath  controller  for  35°  C 

2.  Allow  bath  temperature  to  stabilize. 

3.  Carefully  place  Insert  Board  into  bath  liquid.    EXTREME  care  must  be  taken  to 
ensure  that  none  of  the  wires  are  deformed  or  destroyed. 

4.  Allow  insert  board  to  come  to  thermal  equilibrium  with  the  bath, 

this  takes  approximately  five  minutes. 

5.  Take  temperature  reading  with  platinum  wire  temperature  probe  and  commutator 
bridge,  and  resistance  measurement  of  wires  with  data  acquisition  unit. 

6.  Remove  insert  board,  and  set  bath  temperature  controller  for  a  5°  C  increase. 

7.  Repeat  step  2  thru  6  until  an  90°  C  reading  has  been  achieved. 
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8.  Activate  External  Cooling  system. 

9.  With  temperature  control  setting  decreasing  at  5°  C  intervals,  repeat  step  2  thru  6 
until  a  35°  C  reading  has  been  achieved. 

D.  IN  PLACE  CALIBRATION 

To  ensure  that  the  calibration  curve  of  platinum  wire  1  thru  5  are  reliable,  an  in 
place  calibration  run  was  performed  on  the  wires.  This  calibration  procedure  was  per- 
formed before  any  experimental  runs  were  conducted,  and  after  the  last  experimental  run 
was  completed  to  detect  any  changes  in  calibration  due  to  wire  deformation  or  aging. 

The  equipment  utilized  for  the  inplace  calibration  was  the  same  Data  Acquisition 
equipment  used  for  the  resistance  measurements  of  the  Bath  Calibration  runs  with  two 
additions. 

1.  HP44713A  24  Channel  Multiplexer  with  electronic  cold  junction. 

2.  Four  type  T  thermocouples. 

The  thermocouples  are  now  utilized,  instead  of  the  platinum  resistance  temperature 
probe,  for  measuring  fluid  bulk  temperature.  The  uncertainity  associated  with  the 
thermocouples  is  0.4°C  as  compared  to  the  platinum  resistance  temperature  probe  used 
as  a  standard  Figure  30  shows  the  thermocouples  uncertainity  compared  to  the 
platinum  resistance  temperature  probe. 
1.     Inplace  Calibration  Procedure. 

The  following  procedure  is  used  for  inplace  calibration  of  the  platinum  wires: 

1.  Ensure  Data  Acquisition  system  is  wired  for  4-wire  ohm  measurements. 

2.  Turn  on  Data  Acquisition  system  and  allow  one  hour  of  warm-up  time. 

3.  Load  Calibration  Acquisition  Program  into  the  HP-300  series  Computer.  See  Ap- 
pendix D  for  program  listing. 

4.  Set  power  supply  to  bulk  heaters  to  20volts  and  1.8  amps. 

5.  Start  Acquisition  program.  The  acquisition  program  will  measure  the  resistance 
of  each  platinum  wire,  and  measure  the  temperature  output  of  each  thermocouple 
at  a  10  minute  interval. 

6.  Allow  acquisition  program  to  run  until  bulk  fluid  has  reached  Saturation  temper- 
ature of  56°  C. 

E.  CALIBRATION  RESULTS 
1.     Platinum  Wire  1 

From  Figure  31  the  importance  of  conducting  an  in  place  calibration  prior  to 
any  experimental  runs  becomes  obvious,    wire  l's  Bath  calibration  data  has  been  con- 
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siderably  altered  while  it  was  being  inserted  into  the  test  chamber.  The  initial  in  place 
calibration  data  was  used  for  the  experimental  runs  calibration  curve.  The  resulting 
calibration  was: 

T,  =  53.902  x  R,  -  254.470  (°C) 

From  Figure  32  it  can  be  seen  that  the  post  calibration  curve  matches  well  with 
the  data  acquisition  curve. 

2.  Platinum  Wire  2 

From  Figure  33,  it  is  apparent  that  the  initial  in  place  calibration  curve  fairs  in 
well  with  the  calibration  bath's  decreasing  data.  The  initial  in  place  calibration  data  was 
used  for  the  experimental  runs  calibration  curve.  The  resulting  calibration  was: 

T2  =  53.666  x  R2  -  253.009  (°C) 

From  Figure  34  it  can  be  seen  that  the  Post  run  calibration  matches  well  with 
the  Data  Acquisition  Calibration  curve. 

3.  Platinum  Wire  3 

From  Figure  35  it  is  apparent  that  wire  3's  Bath  calibration  data  has  been  al- 
tered between  the  calibration  bath  and  insertion  into  the  test  chamber.  The  initial  in 
place  calibration  data  was  used  for  the  experimental  runs  calibration  curve.  The  result- 
ing calibration  was: 

T3  =  53.141  xR3-  252.729  (°C) 

From  Figure  36  it  can  be  seen  that  the  Post  run  calibration  matches  well  with 
the  Data  Acquisition  Calibration  curve. 

4.  Platinum  Wire  4 

From  Figure  37,  it  is  apparent  that  the  initial  in  place  calibration  curve  fairs  in 
well  with  the  calibration  bath's  decreasing  data.  The  initial  in  place  calibration  data  was 
used  for  the  experimental  runs  calibration  curve.   The  resulting  calibration  was: 

T4  =  53.484  x  R4  -  252.848  (°C) 

Form  Figure  38  it  can  be  seen  that  the  post  run  calibration  matches  well  with 
the  data  acquisition  calibration  curve. 
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5.     Platinum  Wire  5 

From  Figure  39  it  is  clear  to  see  the  shift  in  calibration  curves  between  the  cal- 
ibration bath  data  and  initial  inplace  calibration  data.  The  calibration  data  has  been 
altered  while  it  was  being  inserted  into  the  test  chamber  obviously.  The  initial  data 
calibration  was  used  for  the  experimental  run  calibration  curve.  The  resulting  cali- 
bration was: 

T5  =  53.085  x  R5  -  260.049  (°C) 

From  Figure  40  it  can  be  further  seen  that  a  considerable  error  has  been  intro- 
duced into  wire  5's  calibration,  an  error  of  3.7°  C.  This  places  an  uncertainity  value  onto 
wire's  5  temperature  data  of  at  least  3.7°  C,  an  unacceptable  value. 
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Figure  30. 


Thermocouple  Uncertainity:  Uncertainity  of  thermocouple's 
tempertaure  measurements  as  compared  to  a  platinum  resistance  tem- 
perature probe. 
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Figure  31.  Calibration  of  Platinum  Wire  1:  Calibration  curve  of  wire  1  of  Cali- 
bration Bath  Temperature  Increasing,  Calibration  Bath  Temperature 
Decreasing,  and  of  the  Initial  in  place  calibration  data. 
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Figure  32.      Calibration  of  Platinum  Wire  1:     Calibration  curve  of  wire  1  oflnitial 
in  place  calibration  curve,  and  of  Final  in  place  calibration  curve 


55 


o 

M 


:  V  i 

z 

w 

5i 

pq 

D 

<! 

z!w 

Sfe 

Mi  . 

i-tH'i— 1 

i 

o!< 
i 

!          :    * 

•  \  \ 

•  \  \ 
:        \ 

:         \ 
; 

:          ;  \  \     : 

\  V 
■           *       v  \ 

j i 

■                               A  \ 

'.                ',                '■  \  \ 
'.               '.                '•     v  \ 
•               •      v  \ 

1 j j....©} 

\ 

i          i 

i         l 

i 

i 

i 

99     r  9     E9     29     T9 


0'9     69     99     L'Q     99     9"9     *'9     E9     29 

(swho)  aoMvisisan 


Q 

S3 


19   09 


Figure  33.  Calibration  of  Platinum  Wire  2:  Calibration  curve  of  wire  2  of  Cali- 
bration Bath  Temperature  Increasing,  Calibration  Bath  Temperature 
Decreasing,  and  of  the  Initial  in  place  calibration  data. 
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Figure  34.      Calibration  of  Platinum  Wire  2:     Calibration  curve  of  wire  2  of  Initial 
in  place  calibration  curve,  and  of  Final  in  place  calibration  curve 
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Figure  35.  Calibration  of  Platinum  Wire  3:  Calibration  curve  of  wire  3  of  Cali- 
bration Bath  Temperature  Increasing,  Calibration  Bath  Temperature 
Decreasing,  and  of  the  Initial  in  place  calibration  data. 
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Figure  36.      Calibration  of  Platinum  Wire  3:     Calibration  curve  of  wire  3  of  Initial 
in  place  calibration  curve,  and  of  Final  in  place  calibration  curve 
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Figure  37.  Calibration  of  Platinum  Wire  4:  Calibration  curve  of  wire  4  of  Cali- 
bration Bath  Temperature  Increasing,  Calibration  Bath  Temperature 
Decreasing,  and  of  the  Initial  in  place  calibration  data. 
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Figure  38.      Calibration  of  Platinum  Wire  4:     Calibration  curve  of  wire  4  of  Initial 
in  place  calibration  curve,  and  of  Final  in  place  calibration  curve 
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Figure 


39.  Calibration  of  Platinum  Wire  5:  Calibration  curve  of  wire  5  of  Cali- 
bration Bath  Temperature  Increasing,  Calibration  Bath  Temperature 
Decreasing,  and  of  the  Initial  in  place  calibration  data. 
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Figure  40.      Calibration  of  Platinum  Wire  5:     Calibration  curve  of  wire  5  of  Initial 
in  place  calibration  curve,  and  of  Final  in  place  calibration  curve 

F.     SUMMARY 

The  following  points  were  brought  out  during  the  calibration  process. 

1.  The  process  of  moving  the  platinum  wires  from  the  calibration  bath  and  installing 
them  into  the  test  chamber  resulted  in  a  shift  in  the  calibration  curves. 
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2.  The  most  probable  cause  of  all  the  calibration  shifts  were  the  lead  wires  in  the  back 
of  the  insert  board  being  pulled  and  adding  an  additional  tension  onto  the  platinum 
wires,  thereby  stretching  them.  The  fact  that  the  calibration  curves  consistently 
shifted  to  an  increase  in  resistance  supports  this  theory. 

3.  The  maximum  difference  observed  between  the  Platinum  Wire  Temperature  Probe 
and  the  Type  T  thermocouples  were  0.4°  C.  Therefore,  a  temperature  uncertainity 
for  the  in  place  calibration  will  be  set  at  0.4°  C. 

4.  The  differences  between  the  wire's  data  acquisition  calibration  curve  and  the  wire's 
fmal  calibration  curve  did  not  exceed  the  thermocouples  uncertainity  in  temper- 
ature. 
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APPENDIX  B.     SAMPLE  CALCULATIONS 

In  these  sample  calculations  the  output  value  of  wire  1  will  be  used. 

A.    DETERMINATION  OF  DERIVED  VALUES 

1.  Determination  of  the  Average  Bulk  Temperature 

The  average  bulk  temperature  of  the  liquid  was  obtained  by  the  arithmetic  av- 
erage of  the  four  immersed  themocouple  outputs. 

_  (T,  +  T2  4-  T3  +  TQ 

*"*  4.0 

(55.835  +  55.884  +  55.988  +  55.932) 

4.0 
=  55.9098°C 

2.  Determination  of  Wire  Current 

The  current  flow  through  the  platinum  wire  is  obtained  by  measuring  the  volt- 
age drop  across  a  precision  2  ohm  resistor. 

.  * Res 

i-PtWire  ~    d 


2Q. 

(0.74605Volts) 
(1.99940hm) 
=  0.2381Ohms 

B.     DETERMINATION  OF  PLATINUM  WIRE  RESISTANCE 

the  resistance  value  of  the  platinum  wires  are  obtained  by  measuring  the  voltage 
drop  across  the  desired  wire  and  dividing  it  by  the  current  flow  through  it. 

D  ^  PtWire 


^PtWlre  ~     i 

lP:Wire 

(1.7233Volts) 
~  (0.2831  Amps) 
=  6.0874  Ohms 

1.     Determination  of  Wire  Surface  Temperature 

The  wire's  surface  temperature  are  obtained  from  the  calibration  curves  shown 
in  APPENDIX  A,  wire  1  will  be  shown  as  an  example. 
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Tptwui  =  53-902  x  ^Ptmre  ~  254.470 
=  53.902  x  6.0874  -  254.470 
=  73.356°C 

2.  Determination  of  Input  Power 

The  input  power  of  the  wires  are  determined  from  the  current  flow  through  the 
wire  times  the  voltage  drop  across  the  wire. 

Powerj  =  Vprwirei  X  ip^ex 

=  1.7233Volts  x  0.2831Ohms 
=  0.48787Watts 

3.  Determination  of  Heat  Flux 

The  heat  flux  dissipated  by  the  platinum  wire  is  determined  from  the  inputted 
power  divided  by  the  wire's  surface  area. 


Hfluxi  = 


Power 


A, 


(0.48787Watts) 

1.5585E-5m2 
=  31,303.44Watts/m2 


where, 


Ax  =  n  x  Dj  x  L, 

=  n  x  5E-5  m  x  0.09922m 
=  1.5585E-5  'm2 

C.     DETERMINATION  OF  FLUID  PROPERTIES 

1.     Film  Temperature 

The  film  temperature  is  the  aritlimetic  average  of  bulk  fluid  temperature,  and  the 
wire's  surface  temperature. 

T  (Tbulk+TptWire) 

1  film  — 


film  2.0 

(54.0°C  +  58.0°C) 

2.0 
=  56.0°C 
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2.     Thermal  Conductivity 

From  Figure  5  of  the  3M  Corporation  Flourinert  Product  Manual  [Ref.    3  :  p. 
15],  the  thermal  conductivity  has  been  determined  to  be: 

(0.6033 -0.001 15  xTmJ 

k  = 


10. 

=  0.05389Watts/m°C 


atT„/m  =  56°C 

3.  Liquid  Density 

Using  the  expressions  in  Table  4B  and  the  constants  presented  in  Table  4C  of 
the  3M  Corporation  Product  Manual,[Ref.  3:  p.  10],  the  liquid  density  has  been  deter- 
mined to  be: 

Pl  =  (1.740  -0.00261  x  1flbn)  x  1000 
=  l,593.8kg/m3 

4.  Kinematic  Viscosity 

Using  Figure  3  in  the  3M  Corporation  Product  Manual,  [Ref.  3:  p.  13]  and  de- 
termining an  expodintial  curve  fit  valid  from  0°C  to  90°C  yields: 

V  =  AeXP(T^  +  2B73.15X) 

A=  1.203952E-8m2/s 
B=  1,058.4109°K 


at  lfllm  =  56°C 


v  =  3.000E  -  7m2/s 


5.     Specific  Heat 

From  Figure  4  of  the  3M  Corporation  Flourinert  Product  Manual,  [Ref.    3 
p.  14],  the  specific  heat  has  been  determined  to  be: 

CP  =  (0.241111  +  3.70374E  -  4  x  Tfllm)  x  4186 
=  1,096.  U/kg°C 

atT^/m  =  56°C 
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6.    Thermal  Expansion  Coefficient 

Using  the  expression  in  Table  4B  and  the  constants  presented  in  Table  4C  of  the 
3M  Corporation  Product  Manual,  [Ref.  3:  p.  10],  the  thermal  expansion  coefficient  has 
been  determined  to  be: 

«  0.00261 


(1.740-0.00261  xTyj/J 
=  0.00164  1/°C 


atT„/m  =  56°C 


7.     Thermal  Diffusivity 


Is  determined  to  be  at  Tfllm  =  56°C 


a  = 


pxC, 


=  3.0848E  -  8m2/s 


8.     Prandlt  Number 


Prandlt  Number  is  determined  to  be  at  Ty,/m  =  56°C 


ft-T 


=  9.725 
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APPENDIX  C.     UNCERTAINTY  ANALYSIS 
A.     UNCERTAINTY  IN  SURFACE  AREA 

coD    =    0.003  mm 

coL    =    0.8  mm 
A     =     rrDL 


where: 

D     = 

Wire  Diameter 

L     = 

Length  of  platinum  wire 

A     = 

Surface  Area  of  platinum  wire 

dA 
3D 

=    (*L) 

dA 

di 


-     (ttD) 


Therefore, 


CO, 


/_a8_V 
V  99.22  J 


0.003  \i 
0.05 


1/2 


CO. 


=    0.061  or  6.1% 


co,     =     9.43529E-7m' 
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B.     UNCERTAINTY  IN  POWER 


a)V2a    =    0.008% +  8/xV     [Ref  13,^.587] 
Mv,ptmre    =    0.008% +  8^V     [Ref  13,/?g.587] 


Mr.iql    ~     ±  0.00 1Q  max  value 


Q      —      V PtWire  x  I 


PtWire 


~       * PtWire  x    ] 


2ii 


2£2 


<3Q 


2ii 


av 


PtWire 


R 


2Q 


_5Q_ 
dR 


PtWire 


R 


2ii 


=     -V 


+ 


2Q, 


2Q. 


PtWire   '         2 
R2Q 


Wo         = 


^2Q  \2    ,    I     VptP 


R 


2Q. 


XVv,PtWlre]    +  1   ~~R X  "^Q  1    +  I  "VP< 


fffre  x    ] 


2fl 


X  CO 


2ii 


2Q. 


1/2 


CO/ 


CO 


V, PtWire    \i 


0) 


V 


+ 


P.20    \2 


+ 


PtWire 


2Q. 


CO 


2Q.     \2 


R 


2Q. 


1/2 


The  worst  case  of  uncertainty  will  occur  when  the  voltages  measured  for  VPtWlrt  and  V2Q 
are  small,  less  than  0.3  Volts. 


cog 

"o" 


(0.008%  +  0.0027%)2  +  (0.008%  +  0.0027%)2  +  (  {  g™°*  {QQ  %  Y 


1/2 


=    0.0524%  or  5.2384E  -  4 
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C.     UNCERTAINTY  IN  HEAT  FLUX 


q   = 


dq 


_Q 

A 


dQ 

is_ 

dA 


<oq    = 


CO, 


J_ 

A 

*2 


i^)2+(iS): 


1/2 


'(*) 


+ 


CO 


<7     \2 


1/2 


=     [(0.061)2  +  (5.2384E-4)2]1/2 

=     6.1002E-2or%  6.1 

D.  UNCERTAINTY  IN  TEMPERATURE 

The  uncertainty  in  the  thermocouple  measurements  (coTC)  is  determined  by  the  un- 
certainty in  the  platinum  wire  temperature  probe  {(oPtProb,)  and  the  error  in  the 
thermocouple's  in  matching  the  probe  temperature  (coA7r). 

q)tc    =     [(coP[Probe)2  +  (coArc)2]1/2 

=    0.400 1°C 

E.  UNCERTAINTY  OF  WIRE  SURFACE  TEMPERATURE 

From  the  calibration  formula 


*>.     -     0.1 -£ 


a)To    =     0.400 1°C 
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^surf    —      z^PtWlre  +  Tq 


=      tVptWire  Tj +  T0  +  <5T 


Where: 

£    =     slope  of  calibration  curve 

T0     =     axis  intercept  for  R  =  0 

Vptwire   =     Voltage  drop  across  Platinum  Wire 

R2ii    =     Resistance  of  in  series  Resistor 

V2ii    =     Voltage  drop  across  series  Resistor 

ST   =     Difference  between  Data  Acquisition  and  Final  Calibration  for  a  particular  data  po 

dTsurf 


8TQ 


8? surf 


dST 


V^surf  R2Q 

=      £ 


^PtWire  V2ii 


V*- surf  VptWire 

=      £ 


^R2ii  V2Q 


<^surf  R2Q 

=      -£Vi 


av2Q    -   ■*'"■»•  vi, 
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0) 


Tsurf 


VptWlre^2Q. 


R 


2ft 


V 


2ft 


«M    +  (1  «ro)    +     CU ^K.P/WrJ    +     £ 


V 


V 


2a 


2ft  \2  2 

V20  / 


1/2 


Using  data  from  Wire  1  at  the  following  point, 


PtWirt 


=     0.71673Volts 


V2CJ    =     0.24529  Volts 

R2n    =     1.9994  Ohms  for  wire  1 

twM    =     53.902°C/Q 

<5T   =    0.02°C 

This  results  in  an  uncertainty  in  temperature  of, 

^Tsurf.wirel     =     0.854°C 
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APPENDIX  D.     DATA  ACQUISITION  PROGRAMS  USED 

A.     PLATINUM  WIRE  RESISTANCE  CALIBRATION  PROGRAM 

The  following  program  was  used  to  obtain  the  Resistance  value's  for  the  platinum 
wire  specimans  used  in  chapter  2.  This  program  is  written  in  HP  basic  language  for  use 
with  controlling  the  HP-3852A  Data  Acquisition  unit. 


10 

DIM  Ohms(0:l 

31  ) 

20 

OUTPUT 

709; 

"INTEGER  I' 

30 

OUTPUT 

709; 

"RST  000" 

40 

OUTPUT 

709; 

"USE,  000" 

50 

OUTPUT 

709; 

"CONF  OHMF' 

60 

OUTPUT 

709; 

"NPLC  1.0" 

70 

OUTPUT 

709; 

"NRDGS  62" 

80 

OUTPUT 

709; 

'DELAY  0. 1 ' 

90 

OUTPUT 

709; 

"RANGE  10" 

100 

OUTPUT 

709 

;  "I  =  99" 

1  10 

OUTPUT 

709 

;  "  I  =  1  +  ] 

[" 

120 

OUTPUT 

709 

;  "MEAS  OHMF.I" 

130 

ENTER 

709;Ohm5(* ) 

140 

Alo-1  . 

E+30 

150 

J  =  1+J 

150 

Ahi-1 . 

E-30 

170 

A=0. 

180 

FOR  K  = 

1  TO 

60 

190 

A=A+0h 

ms<  K 

) 

200 

IF  Ohns(K  )< 

;aio  then 

210 

Alo=0hms(K 

) 

220 

END  IF 

230 

IF  0hms(K  ): 

>Ahi  THEN 

240 

Ahi=0hi 

N5<K 

) 

250 

END  IF 

260 

NEXT  K 

270 

B=A/G0 

, 

280 

PRINT 

"WIRE  SAMPLE  " ; 

J 

290 

PRINT 

"60  RESISTANCE 

AUERAGE  IS  ";B 

300 

PRINT 

"LO  RESISTANCE 

MEASURED  IS  " ;Alo 

310 

PRINT 

"HI  RESISTANCE 

MEASURED  IS  ";Ahi 

320 

IF  J  =  5 

THEI^ 

J 

330 

GOTO  3 

60 

340 

END  IF 

350 

60T0  1 

10 

360 

GOTO  3' 

70 

370 

DIM  Temp5<G 

):61  ) 

380 

K=0. 

390 

OUTPUT 

709; 

"RST  000" 

400 

OUTPUT 

709; 

"USE  000" 

410 

OUTPUT 

709; 

"INTEGER  I 

" 
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420  OUTPUT  709 
430  OUTPUT  709 
440  OUTPUT  709 
450  OUTPUT  709 
460  OUTPUT  709 
470  OUTPUT  709 
480  OUTPUT  709 
490  OUTPUT  709 


"CONF  TEMPT" 
"NPLC  1.0" 
"NRDGS  62" 
"DELAY  0.1" 
"RANGE  AUTO" 
"1=299" 
"1=1+1" 
"MEAS  TEMPT  ,1" 


500  ENTER  709;Temps( * ) 

510  K=K+1 

520  A=0. 

530  FOR  L=1  TO  60 

540  A=A+Temps<L) 

550  NEXT  L 

560  Temp=A/60. 

570  PRINT  "THERMOCOUPLE" ;K 

580  IF  K=4  THEN 

590  GOTO  620 

600  END  IF 

610  GOTO  480 

620  END 


"TEMPERATURE" ; Temp 
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B.     MAIN  DATA  ACQUISITION  PROGRAM 

The  following  program  was  used  to  obtain  the  Temperature  value's  for  the  platinum 
wire  specimans  used  in  chapter  3.  This  program  is  written  in  HP  basic  language  for  use 
with  controlling  the  HP-3852A  Data  Acquisition  unit. 


10 

20 

30 

40 

50 

60 

70 

80 

90 

100 

1  10 

120 

130 

140 

150 

160 

170 

180 

190 

200 

210 

220 

230 

240 

250 

260 

270 

280 

290 

300 

310 

320 

330 

340 

350 

360 

370 

380 

390 

400 

410 

420 

430 

440 

450 


DIM   Te 

L=0 

J  =  0 

L=L+1 

IF    L  =  8 

GOTO    1 

END    IF 

OUTPUT 

OUTPUT 

OUTPUT 

OUTPUT 

OUTPUT 

OUTPUT 

OUTPUT 

OUTPUT 

OUTPUT 

OUTPUT 

OUTPUT 

ENTER 

J  =  J+1 

A=0. 

FOR   K- 

A=A+Te 

NEXT    K 

Tk( J  )  = 

IF   J=7 

GOTO    3 

END    IF 

GOTO    1 

DIM   \)o 

J=0. 

Ct=0. 

OUTPUT 

OUTPUT 

OUTPUT 

OUTPUT 

OUTPUT 

OUTPUT 

OUTPUT 

OUTPUT 

OUTPUT 

OUTPUT 

OUTPUT 

ENTER 

J  =  J  +  1 


mps(0: 16  )  ,Tk( 1  :7  >  ,Tt( 1  : 100)  ,Flux( 1 : 100)  ,Tb( 100) 


1    THEN 
450 


709 
709 
709 
709 
709 
709 
709 
709 
709 
709 
709 


"RST   000" 
"USE    000" 
"INTEGER    I" 
"CONF    TEMPT" 
"NPLC    1.0" 
"NRDGS    17" 
"DELAY    0.1" 
"RANGE    AUTO" 
"1=299" 
"1  =  1  +  1  " 
"MEAS    TEMPT, I 


?09;Tercp5(  *  ) 


1  TO  15 
mps( K  ) 

A/15. 
THEN 
00 

70 

It  5(0:  1  1  )  ,Vk<  1:10) 


709; 
709; 
709; 
709; 
709; 
709; 
709; 
709; 
709; 
709; 
709; 
709;U 


"RST  000" 
"USE  000" 
"INTEGER  I" 
"CONF  DCIT 
"NPLC  1.0" 
"NRDGS  12" 
"DELAY  0.1" 
"RANGE  AUTO' 
"1=99" 
"1=1+1 " 
"MEAS  DCV.r 
olts<*  ) 
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460 

IF    J=5    THEN 

470 

60T0    420 

480 

END    IF 

490 

Ct-Ct+1 

500 

A=0. 

510 

FOR    K-1    TO    10 

520 

A=A+Volts<K  ) 

530 

NEXT   K 

540 

VMCt >=A/10. 

550 

IF    J=6    THEN 

5G0 

60T0    590 

570 

END    IF 

530 

GOTO    420 

590 

DIM   Voltw<0: 1 1  ) 

600 

J-0. 

610 

Ct-5. 

620 

OUTPUT    709';  "RST   000" 

630 

OUTPUT    709;"USE    000" 

640 

OUTPUT    709;" INTEGER    I" 

650 

OUTPUT   709;"CONF    DCU" 

660 

OUTPUT    709; "NPLC    1 .0" 

670 

OUTPUT    709; "NRDGS    12" 

680 

OUTPUT    709; "DELAY    0.1" 

690 

OUTPUT    709; "RANGE    AUTO- 

700 

OUTPUT    709; "1=109" 

710 

OUTPUT    709; "1  =  1  +  1  " 

720 

OUTPUT    709;"MEAS    DCU.I" 

730 

ENTER   709;Voltui(*  ) 

740 

J  =  J  +  1 

750 

IF    J=5   THEN 

760 

GOTO    710 

770 

END    IF 

780 

Ct-Ct+1 

790 

A  =  0. 

800 

FOR   K-1    TO    10 

810 

A=A+Uoltui(K  ) 

820 

NEXT   K 

830 

UMCt  )=A/10. 

840 

IF    J=6    THEN 

850 

GOTO   880 

860 

END    IF 

870 

GOTO   710 

880 

Bt  =  (Tk(  1  HTM  2  HTM  3  HTM  4)  )/4 

890 

V1-VM6  ) 

900 

I1=Uk< 1  )/1 .9994 

910 

V2=VM7> 

920 

I2-VM2  )/2.0076 

930 

U3-VM8) 

940 

I3=Uk(3)/1 .9997 

950 

,V4=VM  9) 

960 

I4-VM  40/2.001  8 

970 

V5-VM  10  ) 

980 

I5»VM5>/2.0025 
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990 

1000 

1010 

1020 

1030 

1040 

1050 

1060 

1070 

1080 

1090 

1  100 

1  1  10 

1  120 

1  130 

1  140 

1  150 

1  160 

1  170 

1  180 

1  190 

1200 

1210 

1220 

1230 

1240 

1250 

1260 

1270 

1280 

1290 

1 


300 
310 


1 


1330 
1340 
1350 
1360 
1370 
1380 
1390 
1400 
1410 
1420 
1430 
1440 
1450 
1460 
1470 
1480 
1490 
1500 
1510 
1520 


P1=I1*U1 

P2=I2*V2 

P3=I3*U3 

P4=I4*U4 

P5=I5*V5 

R1=V1/I1 

R2=V2/I2 

!R3=V3/I3 

|R4=U4/I4 

!R5=V5/I5 

Flux1=P1/1 

Flux2=P2/1 

Flux3=P3/1 

Flux4=P4/1 

Flux5=P5/1 

T1=(53 

T2=(53 

T3=(53 

T4=(53 

T5=(53 


5585284E-5 
5585284E-5 
5710005E-5 
5460563E-5 
5710005E-5 
902123*R1  )-254, 
666456*R2)-253. 
141047*R3)-252. 
483645*R4)-252. 
085508*R5)-2S0. 
Tc  =  (Tk(5)+Tk(6  ))/2. 
PRINT  " 

"  AMBIENT  AIR 
Tk(7)  ,Bt  ,Tc 


46995 

00959 

729131 

848797 

04883 


1  " 

n  « 

"4 
"5 


1 


PRINT 

PRINT 

PRINT  " 

PRINT  "WIRE 

! PRINT  " 

PRINT 

PRINT 

(PRINT 

IPRINT 

(PRINT 

PRINT 

(PRINT 

PRINT 

PRINT 

(PRINT 

(PRINT 

(PRINT 

PRINT  "  " 

Tt<  L  )  =  T1 

Tb(L)=Bt 

Flux(L)=Flux1 

WAIT  100 

GOTO  30 

GOTO  1530 

M-1 

PRINT  " 

IF  M=3  THEN 

GOTO  1530 

END  IF 

PRINT  M,Tt(M) 

M=M+1 

GOTO  1470 

END 


TEMPERATURES" 
BULK  FLUID 


CONDENSOR  SURFACE" 


"  ,L 


TEMPERATURE 
(AMPS) 
,T1  ,11  ,\>\ 
J2.IZ.V2 
"  J3.I3.V3 
"  J4.I4.U4 
"  J5.I5.U5 

HEAT  FLUX 
(WATTS) 
.Fluxl  ,P1  ,R1 
,Flux2,P2,R2 
"  .FLUX3.P3.R3 


CURRENT 
(VOLTS) 


VOLTS  " 
(DEGREE 


C) 


POWER 
(OHMS) 


RESISTANCE 
(WATTS/M2) 


"4"  ,Flux4,P4,R4 
"5"  ,Flux5.P5,R5 


WIRE  TEMP 


,Flux(M)  Jb(M) 


HEAT  FLUX 


BULK  TEMP" 
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C.     NUSSELT  EVALUATION  PROGRAM 

The  following  program  was  used  to  obtain  the  Nusselt  Number  for  the  platinum 
wire  specimans  used  in  chapter  3.  This  program  is  written  in  Fortran  language  for  use 
with  controlling  the  IBM-3803  Main  Frame  Computer. 


FILE:  NUSSELT   FORTRAN   Al 

SUBROUTINE  NUSSELT  (X1,T1,M,Q) 
C 

C   THIS  SUBROUTINE  WILL  CALCULATE  THE  HEAT  TRANSFER  COEFFICIENT  FOR  THE 
C   DATA  OBTAINED  FROM  THE  EXPERIMENTAL  RUNS  BY  USING  THE  CHURCHILL  AND 
C  CHU  CORRELATION 
C 

DIMENSION  Xl(x),Tl(X),Q(X) 

REALX4  X1,T1,Q,K,NU,NUD 

0PEN( 29, FILE='RALY',STATUS=' UNKNOWN ',F0RM=' FORMATTED') 
C 

C   DIAMETER  OF  THE  WIRE  BEING  ANALYSED   IS  NOMINALLY  .05  MM 
C 

D=5.0E-5 
C 

C   GRAVITY  TERM  OF  9.81  M/S2 
C 


G=9.81 
C 

DO  100  1=1, M 
C 

C    FILM  TEMPERATURE  OF  THE  PHYSICAL  PROPERTIES 
C 
CCC    TF  =  (XKI)  +  T1CI)  +  56.0J/2. 

TF  =  (XKI)  +  2.XTl(I))/2. 
C 

C   TEMPERATURE  DIFFERENCE  BETWEEN  THE  WIRE  AND  BULK  FLUID 
C 
CCC   TD  =(X1(I)  -  TKI)  +  56.0) 

TD  =  XKI) 
C 

C      THERMAL    CONDUCTIVITY   OF   FLUID        W/MK 
C 

K    =    (0.6033    -    0.00115*TF)/10. 
C 

C   BULK  MODULUS  OF  THE  FLUID    1/K 
C 

BETA  =  0.00261/(1.740  -  0.00261XTF) 
C 

C   KINEMATIC  VISCOSITY  OF  THE  FLUID    M2/S 
C 

NU  =  1.203952E-8  x  EXP( 1058 . 4109/(TF+273 . 15) ) 
C 

C   SPECIFIC  HEAT  OF  FLUID 
C 

CP  =  (0.24111  +  3.70374E-4  x  TF)x«186.00 
C 

C   FLUID  DENSITY   KG/M3 
C 

RHO  =  (1.740  -  0.00261XTF)X1000.0 
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C   THEMAL  DIFFUSIVITY 
C 

ALFA  =  K/(RHO*CP) 
C 

C   PRANDLT  NUMBER  OF  FLUID 
C 

PR  =  NU/ALFA 
C 

C   RALYEIGH  NUMBER  OF  FLUID   MUST  BE  10xx-5<  RA  <  10XX12 
C 

RA  =  (GXBETAXTDXDXX3. )/<NU*ALFA) 

PRINT  x,  RA 
C 

C    NUSSELT  NUMBER  OF  FLUID       NU=HXD/K 
C 

TERM  =  ((0.518X(RAX*(0.25))x(l+(0.559/PR)XX0.6)xx(-5./12.))XXl5 
|+(0.lxRAXx(l./3.))xxi5.)xx(l./15.) 

NUD=2.0/(LOG(1.+  2. /TERM)) 
C     I 
C 

C   THERORICAL  HEAT  FLUX   WATTS/M2 
C 

Q(I)  =  KXNUDXTD/D 
c100   WRITE(29,500)  TD, RA, NUD, PR, Q(I) 

500   FORMAT  C5CF12. 4,2X) ) 
CLOSE  (29) 
RETURN 
END 
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